Signatures of interband scattering in spectroscopic experiments on MgB 2 
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Within a two-band, strong coupling model we analyze .SiS-tunneling conductance, Raman scat- 
tering and optical conductivity measured on MgBi samples. We find that features observed in 
tunneling and Raman spectroscopy at intermediate energies u = A CT + A T ~ 8 ~ WmeV can be 
consistently explained when incoherent scattering between a and n band is considered. We calcu- 
late the optical conductivity and find that strong coupling effects are important for explaining the 
presently available data, and predict the location and magnitude of interband scattering features 
for this spectroscopic probe as well. 
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PACS numbers: 63.20.Kr, 74.20.-z, 74.25.-q 
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Introduction. The discovery of superconductivity in 
MgB% with a critical temperature Tc = ]jj has 

initiated intense research efforts on elucidating the prop- 
erties of this surprisingly simple compound. The main 
reason behind the relatively high transition temperature 
is still unknown. However, there is a growing consensus 
that strong electron-phonon interaction j| is responsible 
for superconductivity in MgB-i- The large isotope effect 
measured by Bud'ko et al. || for boron (a = 0.26 ±0.03) 
clearly shows that phonons associated with boron vibra- 
tions play a significant role in superconductivity in this 
compound. Tunneling and point contact spectroscopy 
experiments pi [|, 0, fj], ||, ^] reveal a distribution of en- 
ergy gap 2A(0)/Tc between 1.1 and 4.5. Recent experi- 
ments of Giubileo et al. (To) pointed the presence of two 
gaps. Evidence of two gaps were found also in the specific 
heat jllj] and recently in Raman scattering experiments 

& 

According to Liu et al. the four Fermi sur- 

face sheets in MgB2 can be grouped into quasi two- 
dimensional a bands and three dimensional 7r bands, and 
the normal and superconducting properties of MgB2 can 
be described by an effective two-band model. Within 
this model they calculated the coupling constants and 
energy gap ratios in the weak coupling regime. Later 
Golubov et aZ.Q provided the two-band decomposition 
of the superconducting Eliashbcrg functions a 2 Fij(cu) 
(where i and j denote a and tt bands), which describe 
the electron-phonon coupling in MgB2 as function of 
frequency. While electron-phonon interaction was clearly 
investigated in large spectrum of approximations, the role 
of scattering between the bands was not analyzed in de- 
tail. 

We turn our attention to the theory of spectroscopic 
properties of MgBi because, similarly to the situation 
in cuprate superconductors, the most severe constrains 
on theoretical models are imposed by spectroscopic mea- 
surements. In fact, the aim of this paper is to point out 
that a simultaneous explanation of the main features in 
SIS tunnel spectroscopy [|| and Raman scattering Q 
can only be done if one allows for the existence of inco- 



herent (momentum non-conserving) single-particle scat- 
tering events between. We have also analyzed recent op- 
tical conductivity results Q. We find, as shown in detail 
below, that agreement between theory and experiment is 
substantially improved over the weak-coupling analysis 
by considering strong coupling effects alone. We have 
also calculated the effect of incoherent interband scatter- 
ing on such measurements. While the presently available 
data does not allow us to firmly conclude, that there is 
evidence for incoherent interband scattering in the opti- 
cal conductivity data as well, it does not contradict the 
presence of the process either. We nevertheless present 
our predictions for optical conductivity in the presence 
of interband scattering, since we believe that our predic- 
tions for this probe will be useful in analyzing future, 
more accurate measurements. 

Theoretical Model. We will adopt a two-band Eliash- 
berg scheme to calculate the various spectroscopic quan- 
tities. The weak coupling version of the model for two- 
band superconductivity was originally proposed by Suhl, 
Matthias and Walker Jll|. However, due to the lack of 
a genuine two-band superconductor at that time, inter- 
est in this model gradually decreased. More recently, the 
Suhl-Matthias- Walker model was revived and applied to 
cuprates by Kresin[[l6). We have used a similar model Jl7]] 
to calculate the penetration depth in MgBi, and we are 
now applying this framework to derive the spectroscopic 
properties as well. The starting point of our analysis is 
the two-band Hamiltonian H — Hq + Hi + Hi— b- Here 
H.0 is describing the free electrons, 
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In the above expression cJ CT (ck CT ) creates (annihilates) 
electrons in a band and d£ (dko-) are similar opera- 
tors acting in tt band. The last term describes the non- 
interacting phonons with the dispersion given by u/ q . The 
second term in the Hamiltonian, Tij corresponds to the 
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interaction of the conduction electrons with the phonons, 
Hi = Yl ^(k,q)c+ +qCT c kff (6 q + 6i q ) 



k,q,cr 



^(k,q)d + k+qCT d k( r(&q + ^ q ) 



(2) 



k.q.tx 



+ ^(k,q) [c+ +qt7 cM& q + 6± q ) + ft.c.](3) 

k,q,er 

Finally we introduce a new term 7ij_b to account for 
the incoherent (momentum non-conserving) interband 
scattering 



Hi-b = y^(Jk q c kj rf qa - + h.c. 



(4) 



. The effect of this last term is the main focus of our 
study, and we will show below, that it plays an important 
role in explaining the observed features in spectroscopic 
measurements. 

Our approach is based on the self-consistent solution 
of the Eliashberg equations corresponding to the Hamil- 
tonian Self-consistency is reached on the imagi- 

nary axis, and Pade approximation pi] ] is then used to 
analytically continue to the real axis. This procedure 
gives accurately the frequency dependence of the gap 
functions A a (uj) and A-n-fw) and the renormalization fac- 
tors Z a (uS) and Z w (u>). The normal self-energy contribu- 
tion to the i band due to the direct scattering from and 
in the j band can be written as l/(2,Tij)gj(iLu n ) where 
gi(iu) n ) — ijj n j y/w% + Af(uj n ) and the anomalous self- 
energy can be written as X/(2rij)fj(iu) n ) where fi(iui n ) = 
A,;/ y/uJn + Af(uj n ). The structure in A(w) usually gives 
rise to observable variations in the density of states at 
a frequency comparable with the sum of phonon energy 
and gap value. However, in MgBi the phonon mode 
most relevant to superconductivity is Eig — 67meV, and 
observation of density of states features at such high en- 
ergies are difficult to be observed experimentally. For the 
present calculation we model the phonon density with a 
lorentzian: F(u>) = A/((u - uji) 2 + w|) - A/(u>% + wf ) 
for |cj — o>i | < LU3 and zero otherwise. Here the fre- 
quency parameters are: u)\ = 67meV, u>2 = 3.5meV, 
and L03 — SOmeV. With these parameters, and taking a 
cut-off frequency at ui c = 200meV, we fix A — 1.31 by 

imposing J F(tu)cLu = 1. The best fit to the experimen- 


tal temperature dependence of the gap function p2[ and 
with the observed Tc — 39K was found for o? aa = 17.5, 
al„ = 7.0, = 1.0, a\ a = 2.0. The effect of Coulomb 
pscudopotential have been considered through the renor- 
malized parameters /x*, = 0.1 and /i* = for i 7^ j. The 
effect of cross-band interaction can be analyzed in terms 
of Tij = t^ 1 which directly describe the scattering rate. 
This is similar with solving the Eliashberg equation for 
an alloy where similar terms appear due to scattering on 



impurities EfJ. Throughout of our calculations we have 
considered that Ta = (no intra-band scattering) and 
that only IVj 7^ for i ^= i. 

Break Junction Tunneling. We start the comparison 
to experiment by analyzing the break junction SIS tun- 
neling experimental data of Schmidt et al. A care- 
ful analysis of the data reveals the presence a dip at a 
frequency u> ~ A ff + A^ ~ lOmeV and no measurable 
effects are seen near the frequency corresponding to the 
larger gap. This suggests that, due to its 31? nature, only 
the 7r band gives a contribution to the tunneling current 
in break junctions, while the a band, due to its quasi 
2D nature, has little overall contribution. The electron- 
phonon interaction alone does not give a satisfactory ex- 
planation for the above mentioned feature, and a new 
term is needed. This term describes incoherent single 
particle tunneling between bands, and is formally similar 
to the proximity term first proposed by McMillan [ p)] . 
Within this extended Suhl-Matthias- Walker model, the 
superconductivity in 7r band can be thought of as being 
induced by the <r-band. Scattering between bands dimin- 
ishes the larger gap while the smaller one is enhanced. 
More importantly, however, self-energy contributions in 
lowest order in interband scattering will give rise to in- 
termediate states of energies A CT + A T ~ 8 — lOmeV, 
precisely at the voltage where the dip in the tunnel con- 
ductance is seen. 

In order to study the effect of the direct inter-band in- 
teraction we have first calculated the density of states for 
each band by using the well known relation Ng^ (lu) = 

N$(fij)Re ((w - r T ) /y/(u - Tt) 2 - AfR) » where 
Ns(u>) and Nn(u>) are the energy dependent density of 
states of the material in the superconducting and normal 
state. We have calculated first the superconducting den- 
sity of states and then using these results we investigate 
the SIS conductance by calculating first the tunneling 
current by a direct convolution: 
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n F {u + eV)] N { s '\uj)N ( s j) {uj+eV) 
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and then taking the derivative as function of bias. In Eq. 
(||) rip(w) is the Fermi-Dirac function and Ng\u>) with 
i = tt, a is the density of states corresponding to one of 
the bands. The conductance can be found by differentiat- 
ing the tunneling current as function of bias and a direct 
quantitative comparison with the experimental data is 
possible. Fig. |l| represents the results obtained for the 
conductance for the 7r-band together with the experimen- 
tal data of Schmidt et al. Q . The interband self-energy 
leads to modification in the structure of the gap function 
at energies uj — A a + A T . Increasing the scattering rates 
between the bands lead to (i) an enhance of the smaller 
gap corresponding to the tt band and to a decreasing of 
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FIG. 1: Results for normalized 7r-band SIS tunneling con- 
ductance (solid line) compared with the experimental data 
of Schmidt et al. M (open circle symbols). The scattering 
rates used in this calculations are r^ iCT = 2.5 and r CTl7r = 0.5. 
The inset shows the calculated SIS tunneling conductance 
for different scattering rates. 



the er-band gap; (ii) increasing of the dip feature at the 
energy A a + in the SIS conductance spectrum. The 
last effect can be verified experimentally by comparing 
samples of different impurity concentrations. The pres- 
ence of impurities will increase such interband scatter- 
ing processes, and an interesting and somewhat counter- 
intuitive increase in gap size is therefore expected with 
the increase of impurity concentration. This would be a 
qualitative consistency check on our theory. 

Raman Spectroscopy. This is the second spectro- 
scopic experiment that we analyze in the same theoreti- 
cal framework and find an excellent agreement between 
the experiment and our model. Early Raman experi- 
ments performed on MgBi were somewhat inconsistent 
due to differences in sample quality: measurements on 
poly-crystaline samples Jl^] reveal the presence of two 
gaps while single crystal measurements restricted to only 
afe-plane polarized spectra show only the presence of a 
one superconducting gap p3| . Recently a more careful 
analysis resolved this inconsistency and the agreement 
is that the Raman spectrum does reveal the presence of 
two bands We have considered the effect of scatter- 
ing between the bands when calculating the Raman sus- 
ceptibility |25| X(q — * 0,uj). The Raman spectrum was 
calculated separately for each band and then the spec- 
trum was calculated as a weighted sum of contributions 
from each band. The inset of Figj^ presents the normal- 
ized contributions coming from each band. There is a 
complete depletion of the Raman spectrum below 2A, 
for each band and the spectrum has a vertical slope at 
2Aj. The weighting coefficients are 14.61 for 7r-band and 



FIG. 2: Calculated Raman spectra at temperature T = 15K 
(solid line) compared with the experimental data of Chen et 
al. jl^] (open symbols). The dotted line represents the two- 
band BCS like model fit. The scattering rates in these calcu- 
lations were F CTi7r = 0.5 and r^ l<T = 3.0. The inset presents 
the normalized results obtained for the Raman spectrum for 
each band separately. 



4.39 for the cx-band leading to a ratio of 3.32 for the con- 
tributions of each band (in contrast with the 2.97 ratio 
obtained in the Suhl-Matthias- Walker model by Chen et 
al. JT^| . In our calculations the Raman vertex was calcu- 
lated for simplicity in the limit 7(p) = 1. Even in this 
limit the effect of scattering between different bands is 
important, and leads to significant modifications of the 
spectrum. Our results for the Raman spectrum are pre- 
sented in Fig.E and compared to the experimental data 
of Chen et al. |12|. The dashed line is the weak coupling 
BCS two-band model fit proposed by the same authors. 
The significant difference between our results and theirs 
is the presence of the hump at an energy u — A CT + A n , 
a feature that theoretically can be obtained by consid- 
ering only the scattering between the bands, and does 
not appear in the simple Suhl-Matthias- Walker model. 
In our calculations the best fit of the experimental data 
was obtained when r CT 7r = 0.5 and = 3.0. The other 
parameters used in the calculations are r a 2 aa = 17.5, 
= 5.5, a 2 all = 1.0, a^ a — 2.0 and the same values for 
the Coulomb pseudopotential as in the case of analysis 
of the tunneling experiment. 

Optical Conductivity. The last spectroscopic experi- 
ment that we focus on is the optical conductivity, where 
the importance of the scattering between the bands was 
analyzed in the framework of Mattis-Bardeen theory [^6| . 
A quantitative comparison with the experimental data on 
optical conductivity was possible in this case as well, as 
shown in Fig. |^. One immediate conclusion is already 
apparent: the overall fit to the experimental data is much 
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sic effect |27| , and is unlikely to correspond to interband 
scattering. However, we believe that future experiments 
will be able to detect the interband structure we predict 
here. 

Conclusions. We have investigated the role of the in- 
coherent inter-band scattering in the framework of two- 
band strong coupling model. We found results consistent 
with three experimental measurements. The main results 
of our work can therefore be summarized as follows: (i) 
superconductivity in MgB<i can be well explained in the 
framework of the two-band model; (ii) The features ob- 
served in tunneling spectroscopy, and Raman scattering 
at energies ui = A a + A T ~ 8 lOmeV can be consis- 
tently explained by considering the incoherent scattering 
between the bands. While there is no clear feature ob- 
served in optical conductivity at this frequency, we pre- 
dict that such features will eventually be observed, (Hi) 
Overall, the fit to the experimental data is much better if 
a strong coupling theory is used, rather than a two-band, 
weak coupling, constant gap model. ||. 

Acknowledgments. We are especially indebted to Drs. 
H. Schmidt and M. Carnahan for supplying us with tun- 
neling and optical spectroscopy data, respectively. We 
wish to thank them as well as Drs. G. Arnold, G. W. 
Crabtree, K. Gray, M. Iavarone, W.-K. Kwok, G. Kara- 
petrov, A. Koshelev, J. Orenstein, and J. Zasadzinski for 
useful discussions. 



FIG. 3: Real part of conductivity ration a{ '(w)/cr\ '(w) for 
different temperatures (solid lines) from bottom to top T = 
6K, 17.5K, 24K, 27K, 30K compared with the experimental 
data extracted from Ref. [pi (dotted symbols for the same 
temperatures from bottom to top). The dotted line is the 
constant gap model for temperature T=6K. The scattering 
rates are V a ,Tr = 0.5 and rV lt r = 2.5. The inset presents the 
temperature dependence of the gap functions for each band 
separately. 



better in the Eliashberg theory than in the constant gap 
model. Due to the small critical temperature obtained for 
the measured sample Q (T c — 30K) we had to mod- 
ify the phonon coupling constants in order to obtained 
the desired critical temperature. Therefore we consider 
c4 ct = 15.5, a\ v = 5.0, a 2 an = 1.0, a\ a = 2.0, and for 
the scattering rates the best fit with the experimental 
data was obtained for Y a ^ = 0.5 and r OT(CT = 2.5. The 
small feature appearing in the real part of the optical 
conductivity at ui = A CT + A T ~ 8meV (in this case this 
energy is smaller due to the smaller critical temperature 
A„(T = OK) = 2.51 meV, A CT (T = OK) = 5.48 meV), 
is due to the modifications in the frequency dependence 
of the gap functions induced by the effect of inter-band 
scattering in the electronic self-energy. The small feature 
in the data near this frequency appears to be an extrin- 



[1] J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, 

J. Akimitsu, Nature 410, 63, (2001). 
[2] J. Kortus, I. I. Mazin, K. D. Belashchenko, V. P. 

Antropov, and L. L. Boyer, Phys. Rev. Lett. 86, 4656, 

(2001). 

[3] S. L. Bud'ko, G. Lapertot, C. Petrovic, C. E. Cunning- 
ham, N. Anderson, and P. C. Canfield, Phys. Rev. Lett. 

86, 1877, (2001). 

[4] G. Karapetrov, M. Iavarone, W. K. Kwok, G. W. Crab- 
tree, and D. G. Hinks, Phys. Rev. Lett. 86, 4374, (2001). 

[5] G. Rubio-Bollinger, H. Suderow and S. Vieira Phys. Rev. 
Lett. 86, 5582, (2001). 

[6] P. Szabo, P. Samuely, J. Kacmarcik, T. Klein, J. Mar- 
cus, D. Fruchart, S. Miraglia, C. Marcenat and A. G. M. 
Jansen, Phys. Rev. Lett. 87, 137005, (2001). 

[7] A. Sharoni, I. Felner and O. Millo, Phys. Rev. B 63, 
220508, (2001). 

[8] A. Kohen and G. Deuscher, Phys. Rev. B 64, 060506(R), 
(2001). 

[9] H. Schmidt, J. F. Zasadzinski, K. E. Gray, and D. G. 

Hinks Phys. Rev. B 63, 220504(R) (2001). 
[10] F. Giubileo, D. Roditchev, W. Sacks, R. Lamy, D. X. 

Thanh, and J. Klein, S. Miraglia, D. Fruchart, J. Marcus, 

Ph. Monod Phys. Rev. Lett. 87, 177008, (2001). 
[11] F. Bouquet, Y. Wang, R. A. Fisher, D. G. Hinks, J. D. 

Jorgensen, A. Junod and N. E. Phillips, Phys. Rev. Lett. 

87, 047001 (2001). 

[12] X. K. Chen, M. J. Konstantinovic, J. C. Irwin, D. D. 
LAwrie, J. P. Franck, Phys. Rev. Lett. 87, 157002, 



5 



[13] 
[14] 

[15] 

[16] 
[17] 
[18] 



[19] 
[20] 



(2001). 

A. Y. Liu, I. I. Mazin, and J. Kortus, Phys. Rev. Lett. 
87, 08705, (2001). 

A. A. Golubov, J. Kortus, O. V. Dolgov, O. Jepsen, Y. 
Kon g. O. K. Andersen B. J. Gibson, K. Ah, R. K. Kre- 



cond-mat/0111262 



H. Suhl, B. T. Matthias, and L. R. Walker, Phys. Rev. 
Lett. 3, 552, (1959). 

V. Z. Kresin, S. A. Wolf, Phys. Rev. B 46, 6458 (1992). 

C. P. Moca, Phys. Rev. B 65, 132509, (2002). 

R. A. Kaindl, M. A. Carnahan, J. Orenstein, D. S. 
Chemla, H. M. Christen, H.-Y. Zhai, M. Pranthaman, 

D. H. Lowndes, Phys. Rev. Lett. 88, 027003, (2002). 
W. L. McMillan, Phys. Rev. 175, 537, (1968). 

J. P. Carbotte Rev. Mod. P hys. 62, 1027, (199 0); F. 
Marsiglio and J. P. Carbotte |cond-mat/0106143|; D. J. 



Scalapino, J. R. Schrieffer, and J. W. Wilkins, Phys. Rev. 
148, 263, (1966). 
[21] H. J. Vidberg and J. W. Serene J. Low. Temp. Phys. 29, 
179, (1977). 

[22] C. Buzea, T. Yamashita, Supercond. Sci. Tech. 14, R115, 

(2001) and references therein. 
[23] J. W. Quilty, S. Lee, A. Yamamoto, S. Tajima, Phys. 

Rev. Lett. 88, 087001, (2002). 
[24] J. W. Quilty, S . Lee, S. Tajima, A. Yamanaka, 



cond-mat/0206506 



[25] C. Jiang, J. P. Carbotte, Phys. Rev. B 53, 11868, (1996) 

and references therein. 
[26] D. C. Mattis, J. Bardeen, Phys. Rev. Ill, 412, (1958). 
[27] M. Carnahan and J. Orenstein (private communication) 



